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Abstract 

vn . 

CN . The spectrum of the 0"''+, , 1 ^ and 1^" heavy-hght hybrids have been calcu- 

lated in HQET. The interpolated current of the hybrid is chosen as gq'yaG'^^T"'hy(x), 
^ ' gq^ja'JbG'^ T^hv{x) and gqa^aG'^i^T'^hy{x). Some strong couplings and decay widths 

'nI" . of the heavy-light hybrids to B{D)'k are calculated by using the QCD sum rules. The 

mass of O"*""*" hybrid with gluon in TM(1 ) or TE(1^ ) mode is found similar, while 

\^ ' their decay widths are different. A two-point correlation function between the pion and 

^^ ■ vacuum is employed and the leading order of I/Mq expansion is kept in our calculation. 

p^; 1 Introduction 

Oh: 

^ ■ It is almost twenty years to search the exotic hadrons such as the glueballs and hybrids. 

There are some special states which are regarded as the candidates of hybrids, especially the 



p(1400) and p(1600) have been studied widely, but no confirmation has been made so far. 
rS \ Recently, these two special states arouse great interest again. The E852 Collaboration at 

d ' BNL[^ has reported a J'^'^ = 1 ^ isovector resonance p(1405) in the reaction ir^p —>■ rjir^n. 

They also reported the mass 1370 ± IQtfo MeV and width 385 ± 401^^5 MeV. The Crystal 
Barrel Collaboration has also claimed to find an evidence in pp annihilation which may be 
resonance with a mass of 1400 ± 20 ± 20 MeV and a width of 310 ± 50t|J] MeV^. The 
confirmation of these states will provide some evidence for the existence of hybrids. At 
present, all the experiences specialize on the light quark hybrids for some reasons, but it's 
necessary to extend the energy region to hybrids including b or c quark, which is also possible 
in the B or r — C factory. 

Theoretically, the spectrum and decay width of hybrids have been calculated widely with 
many methods such as bag model 0, flux-tube model 0, QCD sum rules 0], lattice and 
other models^. However, there are few works about the spectrum and decay width of 
heavy-light hybrids]^]. HQET has led to much progress in the theoretical understanding of 
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the properties of hadrons[^, one may ask a question whether it is suitable in the heavy-hght 
hybrids. For the heavy-hght hybrids, the sum rules' calculation in full QCD theory in Ref . 
shows that the component of gluon gives a contribution more than l.OGeV to the mass, so 
the "light freedom" seems too heavy to keep the 1/Mq expansion available. The calculation 
of the spectrum of hybrids in HQET will give an answer to the question. Our results show 
that the spectrum of heavy-light hybrids including 6 or c quark are close to those in full 
QCD theory, it is suitable to deal with these hybrids within HQET. Besides, compared with 
bbg and ccg hybrid, the spectrum and decay width of the heavy-light hybrid is easier to be 
calculated in HQET. 

It is interesting for the experimentalists to search the exotic l""*" heavy-light hybrids, so 
theoretical determination of the property of these states is necessary and urgent. In full 
QCD theory, the estimation in Ref. ||^ showed that the sum rule for the mass of 1 '" heavy- 
light hybrids had no platform at all. The masses of 1 '" hybrids were given under the main 
assumption that the contribution of gluon condensate is less than 20% of the bare loop. In 
HQET, the ambiguous situation has been improved greatly. The surface of the A versus 
Borel variable r varies little in a large region, which gives a good platform and determines 
the masses of the hybrids. 

According to the MIT bag model 0, the hybrids with the same J^'^ have different internal 
interactions between the patrons which indicate that they are different states, that is to say, 
the gluon in hybrid can be in different TM(1 ) or TE(1+^) mode. In order to predict the 
properties of them, we should choose suitable generating currents corresponding to these 
states for the calculation. In the case of light quark hybrids[p, these different 0^^ states 
were found to have different masses. In heavy-light hybrids case, the calculation shows that 
the mass split of the O"*""^ heavy-light hybrids with gluon in different mode is not large in the 
1/Mq — y oo limit. The mass of 0"''+ hybrids from two different currents, gq'~faG'^^T"-hy{x) 
with gluon in TM(1 ) mode and gqa^aG'^^T"'hy{x) with gluon in TE(l+~) mode, is found 
similar. 

Though the mass of these two different 0''''^ hybrids is similar, their decay widths to 
B{D)7i final states are found different: the decay widths from current gq'jaG'^^T"-hy{x) are 
about 86 MeV or 16 MeV corresponding to Btt or Dtt final states, respectively, while the 
ones from current gqa^aG'^ T°'hy{x) are 11 MeV or 2.6 MeV. 

The decay widths of the 0"''+ and 1 '" hybrids to B{D)'k have been calculated in Ref. [|T0(| , 
where the decay constants of the hybrids were obtained from the formulae in full theory]^. 
Both the strong couplings and the decay constants are calculated in HQET in this paper, 
and the strong coupling of the 0"'""'" hybrid obtained here is much larger than that in Ref. [^ . 
The large difference is from the values of the decay constant calculated in two ways. 

In the calculation of the strong couplings, two-point function between the pion and vac- 
uum is taken use of insteadly to avoid the ambiguity resulting from the double Borel trans- 
formation for the three-point correlation function and the infrared problem in soft pion limit. 
For convenience, the calculation is kept in the leading order of 1/Mq expansion. 

The paper is organized as follows. The analytic formalism of HQET sum rules for the 
spectrum of hybrid is given in Sec. 2. In Sec. 3, We give the numerical results of the 
spectrum and decay constants of hybrid, the comparison of the spectrum with that in full 
QCD theory is given too. In Sec. 4, with the help of two-point correlation function between 
the pion and vacuum, the analytic formalism of HQET sum rules for the strong couplings 



of hybrids is derived and the numerical resuhs of some decay widths were obtained. In the 
last section, we give the conclusion and discussion. 

2 HQET sum rules for the spectrum of the heavy-hght 
hybrid mesons 

As we know in Ref. [Q, the gluon in hybrid can be in different mode, TM(1 ) or TE(1+^) 
mode. To analyze these different hybrids, we must choose suitable generating current. For 
the spectrum of the 0"*"^ and 1""*" heavy-light hybrids with the gluon in TM(1 ) and 
TE(1"' ) mode, respectively, the interpolated current in HQET is chosen as 

j,{x)=gq-f^G''^^T-K{x). (1) 

where q{x) is the light quark field and hy{x) is the heavy quark effective field, v is the velocity 
of the heavy quark. 

Then, we construct the correlation function as 



where 



n^,H = zydW^^'(0|T{j,(a;),j+(0)}|0) (2) 

LU = 2q ■ V. (3) 



Since the free heavy quark propagator in HQET is / drd^x — vt)—^ and the interaction 

of the heavy quark with the gluon field A^ in the leading order of 1/Mq expansion is ghv-Ah. 
Then under the fixed-point gauge x^A^ = (which will be used throughout this paper), the 
full propagator of the heavy quark {0\T{h{x)h{0)\0) in the leading order of 1/Mq expansion 
is the same as the free one. The freedom of heavy quark can be extracted out of the matrix 
element as a delta function, which facilitates the calculation. 

In the operator product expansion, we keep the perturbative term, two gluon condensate, 
three gluon condensate and two quark condensate. The contribution of mixing condensate 
and higher dimension operators are negligible since their smallness.The Feynman diagrams 
are shown in Fig. 1 , where the double line represents the propagator of the heavy quark. Af- 
ter twice suitable Borel transformation, we obtain the Imlls{uj) and /mHt,(u;) corresponding 
to the scalar and vector contribution, respectively 






^{qq)u^' + ^m{qq)uj'-^{gG'), 
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^{qq)uj' + ^m{qq)uj^-^{gG'), 



where the hght quark mass corrections are considered also in these formulae. 
As for the and l"*"" hybrids, the current was chosen as 

j5>.{x)=gqial5Gl^T-K{x). (5) 

The correlation function expanded is similar to the vector current case except for the opposite 
sign for the contribution of the quark condensates and their spectrum will be determined in 
a similar way. 

For the 0++ hybrid with the gluon in TE(1+^) mode, the following current should be 
used to predict the mass 

j{x)=gqa^^Gl^T-K{x). (6) 

The OPE of the Imll{uj) for this current has been carried out as 

ImUiu) = -^y - m{aG')u + 2aM<iqW - ^{qG^)- (7) 

In the chiral limit, the two gluon and two quark condensates vanish. 

On the phenomenal side, the decay constants of hybrids, Fh± , are defined as 

(0|jgiJ(0+)) = FH+m'Jjlv,, m,\H{l-)) = Fn-mf.e,, (0|j|i7'(0+)) = F^+mJ/,^ (8) 

where mu represent the masses of hybrids, e^ is the polarization vector and the two different 
0++ hybrids with gluon in TM(1 ) or TE(1'^ ) mode is represented as if (0^^) and if'(0"^'^), 
respectively. So the correlation function read 

^ ' (2A-CJ) Le uj'-uo ' 



(2A — u;) Jujc uj' — uj 

where the first term of the right side is the pole term resulting from lowest lying resonance 
contribution and the second term represents the contribution of the continuum state and 
higher resonances, uj^ is the continuum threshold. 

Taking use of the dispersion relations for the correlation function to equate the both sides, 
we obtain 



Fl. 



(2A - uj) 

^ H- 

(2A - uj) 



-i/o-rf.;'^I^2^, (10) 



P'2 

^H+ _ _1 ruJcJ /ImnWl 

■K JO Ul' ~U1 



(2A - uj) 
After the Borel transformation |TT|, they are turned into 

F|±e-2^/^= -l!-^duj'Imn{uj')e~-'l^, (11) 



Table 1: Masses of heavy- light hybrids with different J^^ (GeV). 



jpc 


2Ac 


2Ab 


mniqcg) 


niHiqbg) 


rndfull) 


mh{fuU) 


0++ 


4.4 


4.4 


3.5 


6.9 


4.0 


6.8 


0— 


6.8 


6.8 


4.7 


8.1 


4.5 


7.7 


1-+ 


3.6 


3.6 


3.1 


6.5 


3.2 


6.3 


1+- 


3.8 


3.8 


3.2 


6.6 


3.4 


6.5 


0++ 


4.2 


4.2 


3.4 


6.8 


none 


none 



where T is the Borel transformation variable. So the A can be determined as 



(12) 



After the A has been calculated, the decay constant can be carried out according to (p!0[) . 



3 Numerical results of the spectrum and decay con- 
stants of the hybrids 

In this content, we will give the numerical results of the spectrum and decay constants of 
the hybrids. To proceed the process, the mass of the b and c quark are chosen as 4.7 GeV 
and 1.3 GeV, respectively, the condensates are chosen as 



(0|mgg|0) = -{0.1GeV)\ {0\qq\0) = -{OMGeVf, 

{0\asG^\0) = OMGeV\ {0\gG^\0) = {0.27GeV'^){asG'^). 



(13) 



and the scale of running coupling is set at the Borel parameter T. 

The continuum threshold is chosen as below in the calculation: Uc = 5.0 GeV for the 
and two 0^^ hybrids from current J^(x) and j{x), and Uc = 4.5 GeV for the l"*"" and 1^^ 
hybrids. 

We display our results in HQET and those in full QCD theory in table. 1. In this table, 
the right two columns represent the mass of heavy-light hybrids calculated in full theory, 
the left represent the estimation in HQET and the bottom of this table represents the 0"*""^ 
hybrids with the gluon in TE(1"'"~) mode. 

From the table, the mass of hybrids including b or c quark in HQET is found similar to 
that in full theory, the light freedom in hybrids is not large enough to break down the 1/Mq 



Table 2: Decay constants of heavy-light hybrids ( GeV'^^'^ ). 



hybrid 


Fh+{TM) 


Fh- 


F'h^{TE) 


c quark 


0.31 


0.28 


0.97 


h quark 


0.33 


0.29 


1.01 



expansion. So the calculation in HQET is suitable, which implies that the l/Mg correction 
to the sum rules is not large. 

In 1^+ hybrid case, the sum rule in full theory does not stabilize[^, which may indicate 
that no resonance exists in the channel. The masses of these states were given under the 
main assumption that the gluon condensate contribution is less than 20% of the bare loop. 
In HQET, the 2A of l""*" hybrids versus Borel variable r vary little in a large region, which is 
shown in Fig. 2. The doted line represents that of b quark hybrid and the real line represents 
that of c quark, the little difference between them comes from the running coupling. The 
improvement of the ambiguity may come from the reason|[TT|] that the T, 2A and Uc become 
constant low-energy parameters in the Mq — > oo limit in HQET, while the dependence of 
the parameters M^ and Sc on the heavy quark mass is priori not determined in full theory. 
Besides, the assumption in Ref. 0] is proved reasonable here. 

In the case of light quark hybrids case, the mass of 0"^"^ hybrids with the gluon in different 
mode was found to have a large difference Q. However, the A for the 0"'"+ heavy-light hybrids 
with the gluon in different mode is found similar in the Mq -^ oo limit. The mass of the 
heavy-light hybrids in HQET is represented approximately 



m 



Mq + A + 0(1/Mq), 



(14) 



so the mass split of the 0++ heavy-light hybrids is not large in HQET. The mass of 0+^ 
hybrid with gluon in TM(1 ) mode is about 6.9 GeV and 3.5 GeV corresponding to b or 
c quark hybrid, respectively, and the mass of 0++ hybrid with gluon in TE(1"' ) mode is 
about 6.8 GeV and 3.4 GeV, respectively. 

When the radiative effects are taken account of, the effective current would receive renor- 
malization improvement and the heavy quark expansion of the full current is necessary. 
However, in our derivation, neither the radiative effects nor the 1/Mq correction is taken 
account into. 

The decay constants of the hybrids defined above can be obtained through formula. (p!OD, 
they are all collected in table. 2. The table shows that the decay constants of the 0^+ hybrid 
with gluon in TE(1'^~) mode are larger than those with gluon in TM(1 ) mode. 



Strong couplings and decay widths of heavy-light hy- 
brids 



In Ref. [10|, we have calculated the decay width of 



H,{0+^){k) - 


^BiO-^)ik-q)+n^iq) 


H,{l~^){k) - 


-^5(0-+)(A;-g) + 7r±(g) 



(15) 
(16) 

where the 0"*""^ and 1^+ hybrids with gluon in the TM(1 ) mode and TE(l+~) mode, 
respectively, and the two-point correlation function results from current gq'yaG'^^T°'hy{x) . 
The electric charges of the mesons except for pion have not been written out explicitly. The 
cases of H(.{0^^) -^ D-k^ and Hc{l~^) -^ Dtt"^ have also been calculated there. 

In this section, we will re-consider the same process in HQET firstly. For the decay widths 
of these processes, it is usually calculated through the three-point vertex function or QCD 
light-cone sum rules. However, in order to avoid the ambiguity of the three-point function 
resulted from the double Borel transformation and the infrared divergence in the soft pion 
approximation, we use the following two-point correlator between pion and vacuum 



vv,. 



A,{u',u,v) = iJ dxe"'^7i^{q)\T{ji,{x),j2{0)}\0) = A{u',u)v, + B{u',u){-q, + q 

(17) 
where jiu{x) = gq-y^G1^j^T"'hy{x), J2{x) = hy'j^q^x). A{u',u) and B{u',u) are scalar func- 
tions oi uj = 2k ■ V and uj' = 2{k — q) ■ v , which are determined through the spectral 
density saturated by the mesons corresponding to the interpolated currents, respectively. 
The detailed OPE expansion of the A{uj\uj) and B{uj',uj) has been carried out in Ref. [Il0| . 
For the the 0"''^ hybrid with gluon in TE(l+~) mode, the current jiu{x) in the correlation 
function above should be replaced by j[{x) = gqa^aG'^^T°'hy{x). Then for the processes 
(P!3|), we have another correlation function 

C{u;',u;) = ^ J d'xe''''{n^{q)\T{j[{x),MO)m. (18) 

In the infinite heavy quark mass limit, the following approximate relation 

2A-2A' ^uj-uj' = 2q-v, (19) 

will be used in this paper, where A ~ rriH — Mq and A' ~ rrimeson — Mq. Taking into account 
both the single pole terms and the double pole term in the physical side, we can express 
yl(a;',a;), B{u!',uj) and C{uj',uj) functions of the single variable a;', respectively 

T? -e 1/2 2 

A( i\ — ^ H+ JmgH+m-K^H+'^m . Cq , . 

^^'^>- (2A'-.;0'M3 +2]V3^' ^^"^ 

T? f 1/2 2 

^^'^^= (2A'-^')^Ma +2]V^^' ^^^^ 

rpf £ t 1/2 2 

r^l, ,i\ _ ^ H+ Jm9H+rmT^H+''^m , '^2 ,„„x 



where Cq, Ci and C2 are constants. Fj are decay constants of hybrids defined above, fm are 
decay constants of the B or D meson and gH±m-K refers to strong couphngs, they are defined 
as below 

{<d\jD\D) = -ifDml/M, , (0|jb|5) = -?/Bm|/Mfe, (23) 

(vr±(g)D|i7'(0++)) = g'^.j,^ , (vr±(g)i?|if'(0++)) = ^^,^^, 

(7r±(g)D|i7(0++)) = (7H+D. , (7r±(g)D|i7(l-+)) = gH-D.e ■ q, 

(7r±(g)5|i7(0++)) = gn+B. , {7r^{q)B\H {!-+)) = gns^e ■ q. 

The formulae (^) and (|2^ ) is a little different from those in Ref. [p!0| because of the 
different definition of the decay constants of the hybrids. 

Taking use of the dispersion relation and making Borel transformation on them, we will get 
some equation about the strong couplings. After eliminating the q terms with appropriate 
differentiation, these strong couplings are obtained 

/bf3 

9H^m. = ——^^j^-[2A'A'{r) + Ao]e'^'/\ (24) 

/Vf3 

9H-m. = ^ % ^ [2A'B'ir) + i?o]e^^'/^ 

Fh- Jmmjlf.mf^ 



9'H^m. = - ^, ^ % ^ e'^''^[2A'C{r) + Q 



where Aq and Bq have been given in Ref. ||10||, while C{t) and Cq have the form as 

C{t) = 8y2{3[(mH-m„)62-26i]-(mH-mjFi/r}, (25) 

Co = -8V2(mH-m„)Fi, (26) 



where the parameters hi, Fi have been calculated in Ref.|TO| too. 

Before going on the numerical calculation of the strong couplings and decay widths, it is 
necessary to fixing the parameters firstly. The masses of the heavy quarks and heavy mesons 
have been given in Ref. [0 , the decay constants of B and D mesons are chosen as Ref. [|l^ 



The masses and decay constants of the heavy-light hybrids have been computed above. The 
numerical results of the strong couplings are shown as Fig.3 -^ Fig.5, where the value of 
them is determined around r ~ S.OGeV. They are all displayed in table. 3, where the gH-mn 
is dimensionless. 

To the processes (|15D and (0), the decay widths are given by the following formulae 



T{H{0++) - 


-^m(0"+) +7r) = 


9H,n7v |9| 

Svr mjj 


T{H{1-^) - 


-^m(0"+) + 7r) = 


9l,„^ kP 
247r mjj 



^H '^m 



167r mjj 

[m^ — m^) gHrmr 



(27) 



1927im% 



Then in the Mq — ^ 00, the numerical results of the decay widths read 

T{H{0++) -^ B{0-+) + tt) = SQMeV, T{H{0++) -^ D{0-+) + tt) = IGMeV, (28) 

T{H{l-+)^B{0-+) + 7i) = 2.2MeV, T{H{1-+) ^ D{0-+) + n) = l.OMeV. (29) 

r(ff' (0++) -^ B{0-+) + tt) = llMeV, T{H'{0++) -^ D{0-+) + tt) = 2.6MeV. (30) 



Table 3: Some parameters input and strong couplings of hybrids. ( GeV ). 



hybrid 


Mq 


mm 


Im 


mH(0++) 


m^(l-+) 


m'^(0++) 


gH+rriTT 


dH+m-K 


QH-mn 


b 


4.7 


5.28 


0.18 


6.9 


6.5 


6.6 


8.5 


3.2 


2.8 


c 


1.3 


1.87 


0.19 


3.5 


3.1 


3.2 


2.0 


0.8 


0.8 



Though the decay of hybrids appears to follow the S + P selection rule, which means 
that the decay of hybrids to two S-wave mesons are suppressed fl^, the selection rule is not 
absolute. In the flux tube and constituent glue models, it can be broken by wave function 
and relativistic effects, and the bag model predict that it is also possible that the excited 
quark loses its angular momentum to orbital angular momentum |]l5l, the results obtained 
here support this idea. 

The decay widths of the processes to B{D)'k final states for i/(0^+) with gluon in 
TM(1 ) mode is much larger than those for iJ(l^+) with gluon in TE(1+^) mode, the 
reason is that the final states in the later channels are in the P wave. Besides, the decay 
width of the if (O"*"^) — > Bti obtained here is much larger than that we got in Ref.[]TU|, the 
difference is from the decay constant. The decay constant of the O"*"^ hybrid we got there in 
full theory is much smaller than the one calculated above in HQET. 

Though the mass of these two 0+"^ hybrids is almost the same, the strong couplings 
to pion of them are different, so the decay widths of these two different 0^^ hybrids are 
different. The decay width of the H'{{)^^) —>■ B[D)'k are smaller than those corresponding 
to if (O"*""*") -^ B{D)tt. The physical reason about the difference of decay width of these two 
different O^"*" hybrids is unknown to us yet, but the difference between these states provides 
a nice evidence that the decay property of these two 0"''+ hybrids with the gluon in different 
mode is different. 



5 Conclusion and Discussion 



We calculate the spectrum of the 0+"^, , 1^+ and 1^ heavy-light hybrids with different 
currents in HQET, the results from current gq'yaG'^^T''hy{x) and gq'-yalbGa^T^-h^^x) are 
compatible to those in full QCD theory. The calculation shows that the light freedom in 
heavy-light hybrids is not heavy enough to break down the 1/Mq expansion and it is suitable 
to apply HQET to heavy-light hybrid systems. 

The sum rules for the masses of 1 ^ heavy-light hybrids have no platform at all in 
full theory, so the masses of them were given under some assumptions. The ambiguity of 
these sum rules have been improved in HQET, which suggests the reasonableness of the 
assumptions in Ref.0 in another way. In the calculation, the leading order 1/Mq expansion 
approximation is used and only the first two terms in OPE are kept in our estimate, which 
will bring some errors. Besides, the decay constants will bring in large deviation too, it is 



necessary to determine them more precisely. 

Since the gluon in hybrids can be in different mode, the hybrids with the same J^^ but 
different gluon mode are in fact different states. In the case of light quark hybrids, the two 
different 0+^ states have different masses definitely. In the heavy-light hybrids case, the 
masses of these two different states are found similar in the Mq -^ oo limit, however, the 
calculation shows that the decay widths in the processes of these two hybrids to i?(D)7r final 
states are different. The decay width oi H{0^~^) -^ B{D)tt is found about 86(16) MeV, while 
the decay width of H'{<d++) -^ B{D)tt is only 11(2.6) MeV. 

The strong couplings and decay widths of the heavy-light hybrids in the processes of 
H{1^^) — > B{D)n are calculated too. The large difference of the decay widths calculated 



here from those calculated in Ref.[0 lies on the large difference between the decay constants 
calculated in two different ways. 
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Figure caption 

Fig. 1: Feynman diagrams contributing to the correlation function in HQET. 

Fig. 2: 2A of the l""*" heavy-hght hybrids versus Borel variable T. 

Fig. 3: Strong coupling of H{0^^) heavy-light hybrids versus Borel variable r. 

Fig. 4: Strong coupling of H'{0'^~^) heavy-light hybrids versus Borel variable r. 

Fig. 5: Strong coupling of iJ(l~^) heavy-light hybrids versus Borel variable r. 
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